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^ '. ABSTRACT 

Q ■ Weak gravitational lensing surveys have the potential to directly probe mass density 

fluctuation in the universe. Recent studies have shown that it is possible to model the 
statistics of the convergence field at small angular scales by modeling the statistics 
0^ ' of the underlying density field in the highly nonlinear regime. We propose a new 

O^l i method to model the complete probability distribution function of the convergence 

field as a function of smoothing angle and source redshift. The model relies on a 
hierarchical ansatz for the behavior of higher order correlations of the density field. We 
compare our results with ray tracing simulations and find very good agreement over 
a range of smoothing angles. Whereas the density probability distribution function 
is not sensitive to the cosmological model, the probability distribution function for 
the convergence can be used to constrain both the power spectrum and cosmological 
parameters. 

' Key words: Cosmology: theory - large-scale structure of the Universe - Methods: 

, analytical 
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' Weak distortions in the images of high red-shift galaxies by the intervening large-scale structure provide us with valuable 
^ ] information about the mass distribution in the universe. The study of such distortions provides us a unique way to probe 
• • , the statistical properties of the intervening density distribution. Traditionally, the study of gravitational clustering in the 
. ^ quasi-linear and non-linear regime are done by analyzing galaxy catalogs. However such studies can only provide us with 
, information on how the galaxies are clustered, and to infer the statistics of the underlying mass distribution from galaxy 
j_j ■ catalogs one needs a prescription for galaxy bias. Weak lensing studies are free from such effects and can directly probe the 
■ statistics of the underlying mass distribution. 

Pioneering works in this direction were done by Blandford et al. (1991), Miralda-Escude (1991) and Kaiser (1992) based 
on the earlier work by Gunn (1967). Current progress in weak lensing can broadly be divided into two distinct categories. 
Whereas Villumsen (1996), Stebbins (1996), Bernardeau et al. (1997), Jain & Seljak (1997) and Kaiser (1998) have focussed on 
the linear and quasi-linear regime by assuming a large smoothing angle, several authors have developed a numerical technique 
to simulate weak lensing catalogs. Numerical simulations of weak lensing typically employ N-body simulations, through which 
ray tracing experiments are conducted (Schneider & Weiss 1988; Jarosszn'ski et al. 1990; Lee & Paczyn'ski 1990; Jarosszn'ski 
1991; Babul & Lee 1991; Bartelmann & Schneider 1991, Blandford et al. 1991). Building on the earlier work of Wambsganns 
et al. (1995, 1997, 1998) detailed numerical study of lensing was done by Wambsganns, Cen & Ostriker (1998). Other recent 
studies using ray tracing experiments have been conducted by Premadi, Martel & Matzner (1998), van Waerbeke, Bernardeau 
& Mellier (1998), Bartelmann et al (1998), Couchman, Barber & Thomas (1998) and White & Hu (1999). While a peturbative 
analysis can provide valuable information at large smoothing angles, it can not be used to study lensing at small angular 
scales as the whole perturbative series starts to diverge. 

A complete analysis of weak lensing statistics at small angular scales is not available at present, as we do not have a 
similar analysis for the underlying dark matter distribution. However there are several non-linear ansatze which predicts a tree 
hierarchy for matter correlation functions and are thought to be successful to some degree in modeling results from numerical 
simulations. Most of these ansatze assumes a tree hierarchy for higher order correlation functions and they disagree with each 
other by the way they assign weights to trees of same order but of different topologies (Balian & Schaeffer 1989, Bernardeau 
& Schaeffer 1992; Szapudi & Szalay). The evolution of the two-point correlation functions in all such approximations are left 
arbitrary. However recent studies by several authors (Hamilton et al 1991; Peacock & Dodds 1994; Nityanada & Padmanabhan 
1994; Jain, Mo & White 1995; Peacock & Dodds 1996) have provided us accurate fitting formulae for the evolution of the 
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two-point correlation function, which can be used in combination with these hierarchical ansatze to predict the clustering 
properties of dark matter distribution in the universe. 

Most recent studies in weak lensing have mainly focussed on lower order cumulants (van Waerbeke, Bernardeau k, Mellier 
1998; Schneider et al 1998; Hui 1999; Munshi, Jain & White 1999a), cumulant correlators (Munshi, Jain & White 1999a) 
and errors associated with their measurement from observational data using different filter functions (Reblinsky et al. 1999) . 
However it is well known that higher order moments are increasingly sensitive to the tail of the distribution function and are 
more sensitive to measurement errors due to the finite size of a catalog. On the other hand numerical simulations involving 
ray tracing experiments have already shown that while the probability distribution function associated with the density field 
is not sensitive to cosmological parameters, its weak lensing counterpart can help us in the estimation of such parameters 
from observational data (Jain, Scljak & White 1999). At present there is no prescription for the theoretical estimation of 
the PDF for the smoothed convergence field «;(^o)- Valageas (1999) has used a hierarchical ansatz for computing the error 
involved in the estimation of f2o and A from SNela observations. A similar fitting function was recently proposed by Wang 
(1999). Our formalism is similar to that of Valageas (1999) although the results obtained by us include the effect of smoothing. 
The formalism we introduce can be easily extended to the multi-point PDF and hence to compute bias and Sn parameters 
associated with "hot spots" of convergence maps. Analytical results and detailed comparison against numerical simulations 
will be presented elsewhere. 

The paper is organized as follows. In section 2 we briefly the ray tracing simulations. Section 3 presents the analytical 

results necessary to compute the PDF of the smoothed convergence field k(6q). The comparison of these results with ray 
tracing simulations is made in section 4. Section 5 contains a discussion of our results. 



2 GENERATION OF CONVERGENCE MAPS FROM N-BODY SIMULATIONS 

Convergence maps arc generated by solving the geodesic equations for the propagation of light rays through N-body simulations 
of dark matter clustering. The simulations used for our study arc adaptive P'^M simulations with 256'* particles and were 
carried out using codes kindly made available by the VIRGO consortium. These simulations can resolve structures larger than 
30h~^kpc at « = accurately. These simulations were carried out using 128 or 256 processors on CRAY T3D machines at 
Edinburgh Parallel Computer Center and at the Garching Computer Center of the Max-Planck Society. These simulations 
were intended primarily for studies of the formation and clustering of galaxies (Kauffmann et al 1999a, 1999b; Diaferio et 
al 1999) but were made available by these authors and by the Virgo Consortium for this and earlier studies of gravitational 
lensing (Jain, Seljak & White 1999, Reblinsky et al. 1999, Munshi & Jain 1999a) 

Ray tracing simulations were carried out by Jain et al. (1999) using a multiple lens-plane calculation which implements 
the discrete version of recursion relations for mapping the photon position and the Jacobian matrix (Schneider & Weiss 1988; 
Schneider, Ehlcr & Falco 1992). In a typical experiment 4 x 10*' rays are used to trace the underlying mass distribution. The 
dark matter distribution between the source and the observer is projected onto 20 - 30 planes. The particle positions on each 
plane are interpolated onto a 2048^ grid. On each plane the shear matrix is computed on this grid from the projected density 
by using Fourier space relations between the two. The photons are propagated starting from a rectangular grid on the first lens 
plane. The regular grid of photon position gets distorted along the line of sight. To ensure that all photons reach the observer, 
the ray tracing experiments are generally done backward in time from the observer to the source plane at red-shift z = Zg. 
The resolution of the convergence maps is limited by both the resolution scale associated with numerical simulations and also 
due to the finite resolution of the grid to propagate photons. The outcome of these simulations are shear and convergence 
maps on a two dimensional grid. Depending on the background cosmology the two dimensional box represents a few degree 
scale patch on the sky. For more details on the generation of /t-maps, see Jain et al (1999). 



3 ANALYTICAL PREDICTIONS 

In this section we provide the necessary theoretical background for the analytical derivation of the probability distribution 
function for the smoothed convergence field k(^o)- We will use models of nonlinear hierarchical clustering; details of various 
hierarchical ansatzes can be found in Balian & Schacffcr (1989), Bernardeau & Scliacffcr (1992,1999), Bernardeau (1992,1994), 
Szapudi & Szalay (1993), Munshi et al. (1999a,b,c), Valageas & Schaeffer (1998), Valageas et al. (1997)). 
The following line element describes the background geometry of the universe: 

dr^ = -cdf + a{t){dx' + r^(x)d^f^) (1) 

Where we have denoted the angular diameter distance by r(x) and scale factor of the universe by a{t)\ r(x) = K~^^^ 8m.{K~^^'^x) 
for positive curvature, r(x) = (— Jf)"^^^ sinh((— i<r)~^^^x) for negative curvature and r(x) = x for a zero curvature universe. 
The curvature K is given in terms of the present value Hq and V,o as K = (fio — i)H§. Parameters charecterising different 
cosmological models that we will be studying are listed in Table 1. 



Table 1. Cosmological parameters characterizing diflferent models 
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3.1 Formalism 

The statistics of the convergence k in weak lensing is very similar to that of the projected density. In our analysis we will 

consider a small patch of the sky where we can use the plane parallel approximation or small angle approximation to replace 
spherical harmonics by Fourier modes. The 3-dimensional density contrast S along the line of sight when projected onto the 
sky with the weight function a;(x) gives the convergence in a direction 7. 

«(7) = / dxuj{x)Sir{xh) (2) 
Jo 

In our discussion we will be assume the source gala^Kies to be at a fixed red-shift Zg. For this case, one can express the weight 
function as, 

^(X)-^— , (3) 

where Xs is the comoving radial distance to the source galEixies. Using the Fourier transform of d, the convergence can be 
expressed as: 



«:(7) = / dxi^ix) / ——jexp{ixk\i+irek±)5k, (4) 
Jo J (^tt) 

where 6 denotes the angle between the line of sight direction 7 and the wave vector k, fey and k± denote the components of 
k, parallel and perpendicular to the line of sight. In the small angle approximation one assumes that A;x >> k\\. Using these 
definitions we can compute the projected two-point correlation function in terms of the dark matter power spectrum P{k, x) 
(Peebles 1980, Kaiser 1992, Kaiser 1998): 



i^M^M) = ['dX^ I exp(W) H^yX). (5) 

Where we have introduced 1 = r(x)k_L to denote the scaled wave vector projected on the sky. The variance of k smoothed 
over an angle 9o is (Jain & Seljak 1997) 

^^(^0)) = jyx.-^ J ^ ^{^yx)wiiieo) (6) 

Similarly the higher order moments of the smoothed convergence field relate {k''{9o} to the 3-dimensional multi-spectra of the 
underlying dark matter distribution Bp (Hui 1999, Munshi & Coles 1999a): 



= f .X.^/|i^M..(..o)/|i^H^.(..o)/0^^ (7) 

(«^(^o)) = _^"dxi^ I ^w.iMo) I ^w.iMo) I ^w.iiseo) I ^w.iueo) 



E' 

(8) 
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'V^'Kx)'^'^'^Vi,=o- 



We will use these results to show that it is possible to compute the complete probability distribution function of k from 
the underlying dark matter probability distribution function. Details of the analytical results presented here can be found in 
Munshi & Coles (1999b). 



3.2 Hierarchical Ansatze 

The spatial length scales corresponding to small angles are in the highly non-linear regime. Assuming a tree model for the 
matter correlation hierarchy in the highly non-linear regime, one can write the general form of the iVth order correlation 
function as (Groth & Peebles 1977, Pry & Peebles 1978, Davis & Peebles 1983, Bernardeau & Schaeffer 1992, Szapudi & 
Szalay 1993): 



(JV~1) 

U(ri,...rN)= ^^-^ n ^('■■''■J)- (9) 

a, N — trees labcllings edges 

It is interesting to note that a similar hierarchy develops in the quasi-linear regime in the limit of vanishing variance 
(Bernardeau 1992); however the hierarchal amplitudes Qn.c become shape dependent functions in the quasilinear regime. 
In the highly nonlinear regime there are some indications that these functions become independent of shape, as suggested 
by studies of the lowest order parameter Q3 — Q using high resolution numerical simulations (Sccocimarro et al. 1998). In 
Fourier space such an ansatz means that the hierarchy of multi-spectra can be written as sums of products of the matter 
power-spectrum. 

B2(kl,k2,k3)^^^^^0 = 0(P(kl)P(k2) +P(k2)P(k3) +P(k3)P(kl)) (10) 

B3(ki, ka, ka, k4)^ = Ra P(ki)P(ki + k2)P(ki + k2 + k3) + cyc.perm. + Rt P(ki)P(k2)P(k3) -I- cyc.perm. 



Different hierarchical models differ in the way they predict the amplitudes of different tree topologies. Bernardeau & Schaeffer 
(1992) considered the case where amplitudes in general are factorizable, at each order one has a new "star" amplitude and 
higher order "snake" and "hybrid" amplitudes can be constructed from lower order "star" amplitudes (see Munshi, Melott & 
Coles 1999a, b,c for a detailed description). In models proposed by Szapudi & Szalay (1993) it was assumed that all hierarchal 
amplitudes of any given order are degenerate. Galaxy surveys have been used to study these ansatze. Our goal here is to 
show that weak-lensing surveys can also provide valuable information in this direction, in addition to constraining the matter 
power-spectra and background geometry of the universe. We will use the model proposed by Bernardeau & Schaeffer (1992) 
and its generalization to the quasi-linear regime by Bernardeau (1992, 1994) to construct the PDF of the convergence field 
k{6o). We express the one-point cumulants as: 



{^'(7)) = {3Q3)C3[4o] = Ss{K'(eo)f (11) 

{«;*(7)) = {12Ra + ARt)C44o] = S4{ti\eo)}\ (12) 
where 

C4-'"(eo)] = ^^0^^nZdx, (13) 
and 

^oo = / lii^P (-I-t) WiiWo). (14) 



(27r)2 yrix) 

Equation ( |ll| ) was derived by Hui (1998). He showed that his result agrees well with the ray tracing simulations of Jain, Seljak 
and White (1998). More recent studies have shown that two-point statistics such as cumulant correlators can also be modeled 
in a similar way (Munshi & Coles 1999, Munshi, Jain & White 1999). 



3.3 Scaling Models and Generating Functions 

The success of analytical results on lower order cumulants, when compared with numerical ray tracing experiments, motivates 
a more general analysis of the probability distribution function of the smoothed convergence field k{6o). For this purpose we 
found the formalism developed by Balian & Schaeffer (1989) and Bernardeau & Schaeffer (1992) to be most suitable. These 
results are based on a general tree hierarchy of higher order correlation functions and the assumption that the amplitudes 
associated with different tree-topologies are constant in the highly non-linear regime. These results were generalized by 
Bernardeau (1992, 1994) to the quasi-linear regime where the perturbative dynamics can be used to make more concrete 
predictions. In this section we review the basic results from the scaling models before extending such models to the statistics 
of the smoothed convergence field n{9o) (for more details see Munshi & Coles 1999b). We will use the small angle approximation 
throughout our derivation. Our results can be generalized to the case of the projected clustering of galaxies. 

In a scaling analysis of the probability distribution function (PDF), the void probability distribution function (VPF) 
plays a fundamental role. The VPF can be related to the generating function of the cumulants or Sn parameters, (j}{y) (White 
1979, Bahan & Schaeffer 1989). 

P40) = exp(7Va(7Ve)) = exp ( - , (15) 

where P„(0) is the probability of having no "particles" in a cell of of volume v, N is the average occupancy of the particles 
and Nc = N£^2- The VPF is meaningful only for a discrete distribution of particles and can not be defined for smoothed 
fields such as S or k{6o). However the scaling functions defined above as u{y) — —(j>{jj)/y, are very useful even for continuous 
distributions where they can be used as a generating function of the one-point cumulants or the Sn parameters. 



p=l 



The function <j)iy) satisfies the constraint 5*1 = 52 = 1 necessary for tfic normalization of the PDF. The other generating 
function which plays a very important role in such an analysis is the generating function for vertex amplitudes Vn, associated 
with nodes appearing in the tree representation of the higher order correlation hierarchy (Q3 = 1/2, Ra = I'i and Rb = 1^3)- 

e(r) = l-r+|.^-|.^ + ... (17) 

A more specific model for Q{t) is sometimes used, which is useful in making more detailed predictions (Bernardeau & Schaeffer 
1979): 

e(r)= (l + f )"'^" (18) 

We will relate k with other parameters of scaling models. While the definition of VPF does not use any specific form of the 
hierarchical ansatz, writing the tree amplitudes in terms of the weights associated with the nodes is only possible when one 

assumes a factorizablc model for the tree hierarchy (Bernardeau & Schaeffer 1992). Moreover, in such an ansatz the generating 
functions for the tree nodes can be related to the VPF by solving a pair of implicit equations (Balian & Schaeffer 1989). 

Hy) = yQ{r) ^ \yr^Q{r) (19) 

T = -y^Q{r) (20) 
The VPF and PDF can be related to each other by the following equation (Balian & Schaeffer 1989) 

P(5) = r ^ exp r (i±%^^l (21) 

So it is clear that the properties of 4>{y) completely determine the behavior of P{S) for all values of 6. However different 
asymptotic expressions of 4>{y) govern the behaviour of P{5) for different intervals of 5. For large y we can express 4>{y) as 

<P{y) = ay'-". (22) 

We have introduced a new variable u) for the description of VPF. This parameter plays a very important role in the scaling 
analysis. No theoretical analysis has been done so far to link u) with the initial power spectral index n. Numerical simulations 

are generally used to fix cj for different initial conditions. Such studies have confirmed that an increase of power on smaller 
scales increases the value of u). Typically an initial power spectrum with spectral index n — —2 (which should model CDM-like 
spectra that we have considered in our simulations at small length scales) produces a value of w = 0.3 which we will use in 
our analysis of the PDF of the convergence field k(0o)- 4>{y) exhibits a singularity for small but negative value of ys 

Hy) = <t>s - asT{ujs){y - ysy^" ■ (23) 

The parameter ka which we have introduced in the definition of G{t) can be related to the parameters a and lo appearing 
in the asymptotic expressions of 4>{y) (Balian & Schaeffer 1989. Bernardeau & Schaeffer 1992). 



W = ka/{ka + 2) (24) 

ka + 2^fc./fc„+2 ^25) 



2 

Similarly the parameters ys and as which describe the behavior of <j}{y) near the singularity can be related to the behavior of 
G(t) near its singularity Ta, which is given by (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992) 

The parameter ys can be expressed as (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992) 
which in turn can be related to ka as 

1 uk.+2r^+' 

ys~ k^{ka + l)'^^+'- ^ ' 

These singularities describe P{5) for large values of 6. The newly introduced variable Xi, will be useful later to define the 
large 5 tail of P{5). The different asymptotes in (f>{y) are linked to the behavior of P{S) for various values of 5. However for 
very large values of the variance i.e. I2 it is possible to define a scaling function P{S) = h{x)/Q which will encode the scaling 




Figure 1. Analytical predictions for the variance of the reduced convergence field t; as a function of red-shift. Note that whereas the 
variance of k{8o) increases with increasing source red-shift, the trend is reversed for r]{6o). Both r]{9o) and k{9o) decrease with increasing 
smoothing radius Oq. 



behavior of the PDF, where x plays the role of the scaling variable and is defined as (1 + 5)/£,2- We list the different ranges 
of 5 and specify the behavior of P{5) in these regimes (Balian & Schaeffer 1989). 



<< 1 + 5 << ^2 



1 + 5 » 6 



P{5) 



'1 + 5 



a I — uj f I + S 



1+5^ 



^fi+i)expf-ii^) 



(29) 



(30) 



The integral constraints satisfied by the scaling function are Si = xh{x)dx — 1 and S2 ~ x'^h{x)dx — 1. These 
constraints lead to the correct normalization of P{&)- Several numerical studies have been conducted to study the behavior 
of h(x) for different initial conditions (e.g. Colombi et al. 1996; Munshi et al. 99). For very small values of 5 the behavior of 
P(S) is determined by the asymptotic behavior of </>(?/) for large values of y. Hence the PDF can be expressed as (Balian & 
Schaeffer 1989): 



P{5) 



dy 
2ni 



exp 



(1 + 5)y - ay~ 
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(31) 



The above equation implies that it is possible to define another scaling function g[z), completely determined by uj, with the 
scaling parameter z = (l + (5)a" 

compared to that from g{z) (e.g. Bouchet & Hernquist 1992). 



-1/(1 However, numerically it is much easier to determine to from the study of (j{y) 



l + 5<<6; P{5)^a~ 



— exp{zt - t ) 



P{S)^a^f^ 



27ra;z(i+'^)/'^ 



exp 



(32) 
(33) 



(34) 



To summarize, the entire behavior of P{5) is encoded in the two scaling functions, h{x) and g{x). These functions are relevant 
for the behavior of P{5) at small and large 5, respectively. Typically P{5) shows a cutoff at both large and small values of 5 
and exhibits a power-law in the middle. The power law behavior is observed when both g{z) and h{x) overlap, in the highly 




Figure 2. The variance of the convergence field (ft2(eo)) as a function of red-shift Zs ■ Difi^erent panels correspond to different smoothing 
angles. Note that the ordering of the variance for the convergence field of different cosmological models is not the same as that of the 
reduced convergence field r}. 



non-linear regime and with the decrease in ^2 the range of 5 for which P{5) shows such a power law behavior decreases, finally 
to vanish in case of small variance i.e. in the quasi-linear regime. 

In the quasi-linear regime exactly the same formalism can be used to study the PDF. However the generating function now 
can be explicitly evaluated by using tree- level perturbative dynamics (Bernardeau 1992). It is also possible to take smoothing 
corrections into account in which case one can have an explicit relation of uj and the initial power spectral index n. In general 
the parameters ka or uj charecterising the VPF or the CPDF are different from their highly non-linear counterparts. 

The PDF now can be expressed in terms of Gs{t) (Bernardeau 1992, Bernardeau 1994): 



P(S) dS = 
Gs{r 



1/2 



exp 



V 2^2 



dr 



(35) 
(36) 



The above expression for P{S) is valid for 5 < Sc where Sc is the value of S which cancels the numerator of the prefactor of 
the exponential function appearing in the above expression. For 5 > 5c the PDF develops an exponential tail which is related 
to the presence of a singularity in 0(y), as was the case of its highly non-linear counterpart (Bernardeau 1992, Bernardeau 
1994), 



P{5) dS - 



-5/2 



exp 



'6 6 



dS 



(37) 



Since the value of the parameter u is different in the two regimes, the associated parameters such as ys and as will also be 
different in these two regimes. We will show that the expressions for the probability distribution function P(ry) of the reduced 
convergence rj will have exactly the same asymptotes as the weakly non-linear PDF of P(l -I- 5), although the parameters 
associated with P(r?) will correspond to the ones used in the highly non-linear regime for P{l + 5). This is one of the important 
result of our analysis. 

It may be noted that similar expressions can be derived for the approximate dynamics sometimes used to simulate 
gravitational clustering in the weakly non-linear regime, e.g. Lagrangian perturbation theory which is an extension of the 
Zeldovich approximation (Munshi et al. 1994). 
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Figure 3. The minimum value of k(6o), i.e. —Kmin is plotted as a function of the source red-shift Zs. The minimum value of k occurs in 
those line of sight direction where smoothed density perturbation (5 = —1 all along the line of sight. It is independent of the smoothing 
angle Bq and depends only on the background geometry. 



3.4 PDF of the Smoothed Convergence Field 

For computing the probability distribution function of the smoothed convergence field n{Oo), we will begin by constructing 
its associated cumulant generating function ^i+K(eo){y)'- 



(38) 



Now using the expressions for the higher moments of the convergence in terms of the matter power spectrum (equations ^ 
and|l|) 
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We can now use the definition of cj>(y) for the matter cumulants to express $i+k(£i„) (y), in terms of 4'{y): 

{r\x)4) 1 J-(x)/(^P(KB)W^'W- 



Ay) 



dx 



r2(x) 



y I ^{x)dx 



(39) 



(40) 



The extra term comes from the A*' = 1 term in the expansion of $1+^(8^). Note that we have used the fully non-linear 
generating function cf) for the cumulants, though we will use it to construct a generating function in the quasi-linear regime. 
The analysis becomes much easier if we define a new reduced convergence field ri{6o): 



= 1 + 



(41) 



where the minimum value of k{6o) i.e. Kmin occurs when the line of sight goes through regions that are completely empty of 
matter (i.e. S — —1 all along the line of sight): 

rxs 

Kmin = dx^(x)- (42) 

Jo 

Figures 1 and 2 compare the variance of rj and k. While k{9o) depends on the smoothing angle, its minimum value 
Kmin depcuds ouly on the source red-shift and background geometry of the universe, independent of the smoothing radius. 
Figure 3 shows the dependence of Kmin on source red-shift for three cosmological models. With the reduced convergence rj, 
the cumulant generating function is given by. 




Figure 4. Analytical predictions for the probability distribution function as a function of the reduced convergence rj. Curves from 

left to right correspond to decreasing values of the variance y'^ = .25, .5, and 1.. In all our comparison with ray tracing experiments 
we have used lo = 0.3. However there are some uncertainities in estimation of uj from numerical simulations (see Colombi et. al. 1996; 
Munshi et al. 1999). In the right panel we show how small changes in u! change the probability distribution function for the reduced 
convergence 7]. 
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The constructed cumulant generating function $17(2/) satisfies the normalization constraints Si = ^2 = 1. The scaling 
function associated with P{ri) can now be easily related with the matter scaling function h{x) introduced earlier: 
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Using this definition we can write: 



(44) 
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Kx)dx 



Luix) J d^k^P{k^)WHk^eor{x)) Jo' dxcoix) 
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Lo(x) J d2kxP(k)M/2(kxeor(x)) J^' rfx^(x) 



While the expressions derived above are exact and are derived for the most general case using only the small angle approxima- 
tion, they can be simplified considerably using further approximations. In the following we will assume that the contribution 
to the X integrals can be replaced by an average value coming from the maximum of (^{x)y i-e. Xc (0 < Xc < Xs)- So we replace 
/ fix)dx by l/2f{xc)^x where is the interval of integration, and /(x) is the function of x under consideration. Similarly 
we replace the uj{x) dependence in the k integrals by cj(xc). Under these approximations we can write, 



-^viv) = (46) 
h,,{x) = h{x) (47) 

Thus we find that the statistics of the underlying field 1 + 5 and the statistics of the reduced convergence -q are exactly the 
same under such an approximation (the approximate functions and hn{x) do satisfy the proper normalization constraints). 
Although it is possible to integrate the exact expressions of the scaling functions, there is some uncertainty involved in the 
actual determination of these functions and associated parameters such as to, ka, from N-body simulations (e.g. see Munshi 
et al. 1999, Valageas et al. 1999 and Colombi et al. 1996 for a detailed description of the effect of the finite volume correction 
involved in their estimation). We have used $17(2/) as derived above to compute Pirj) with the help of equation (pl|). 

An Edgeworth expansion of the pdf can be made by starting with the simpler form P{ri{9o)), which can then be used to 
construct the Edgeworth series for P{k{9o)). The Edgeworth expansion (see e.g. Bernardeau & KofTman 1994) is meaningful 
when the variance is less than unity, which guarantees a convergent series expansion in terms of Hermite polynomials Hniy), 




Figure 5. Analytical predictions for probability distribution function of the weak lensing convergence field P(re) smoothed with an angle 
00 = 1' is plotted as a function of k. The source red-shift is fixed at Zg = 1. The four curves correspond to the SCDM, LCDM, OCDM 
and TCDM models as indicated. The Solid lines with black dots are the results from ray tracing simulations. 



of order n and with u = v/yi&i)- Hence it is used in quasi- linear analysis with the perturbative expressions for cumulants. 
However the parameters used in the expansion of P{ri) are from the highly non-linear regime; i.e. although the variance is 
smaller than unity, the parameters that characterize it are from the highly non-linear dynamics of the underlying dark matter 
distribution. 



1 



(48) 



The magnification ^ can also be used instead of k using the weak lensing relation ^i{9o) — 1 + 2k(6'o). Its minimum value 
can be related to ftmin defined earlier as fj,min = 1 + 2Kmin ■ Finally, the reduced convergence i] and the magnification /i can 
be related by the following equation (Valageas 1999) : 

Vieo) = ^^f°2"^""" . (49) 

So we can express the relations connecting the probability distribution function for the smoothed convergence statistics k{6o), 
the reduced convergence ri{do) and the magnification /i(6'o) as, 

p(K{eo)) = 2P{f,)^PM—^ — -^p{^)-L^. (50) 

Throughout our analysis we have used a top-hat filter for smoothing the convergence field, but our study can be extended 
to compensated filters (Schneider et al. 1998; Reblinsky et al. 1999), which may be more suitable for observational purposes. 
The formalism which we have developed for one-point statistics such as the PDF and the VPF can also be extended to 
compute the bias and higher order cumulants associated with spots in k, maps above a certain threshold. The statistics of 
such spots can be associated with the statistics of over-dense regions in the underlying mass distribution which represent the 
collapsed objects. A detailed analysis of these issues will be presented elsewhere (Munshi & Coles 1999b). 



4 COMPARISON WITH RAY TRACING SIMULATIONS 

To compare analytical results with numerical simulations we smooth the convergence maps generated from numerical simula- 
tions using a top-hat filter of suitable smoothing angle Oq. The minimum smoothing radius we have used is 1' which is much 
larger than the numerical resolution length scale. The maximum smoothing radius we have studied is 8' which is much smaller 




Figure 6. As in Figure 5, with smoothing angle Bq = 2'. 




Figure 7. As in Figure 5, with 



smoothing angle Bq =4'. 
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Figure 8. As in Figure 5, with smoothing angle = 

than the box size and we expect that finite volume effects are not significant in our studies. The box size is L = 166.28' for 
the EDS models, L = 235.68' for the Q. = 0.3 open model, and L = 209.44' for the Q. = 0.3 model with cosmological constant 
A = 0.7. Numerical outputs from one realization of each of these cosmological models were used to compare against theoretical 
results. Elsewhere, we have used these simulations to test analytical results for cumulants and cumulant correlators and found 
very good agreement. 

Figures 4-8 show the analytical and numerical pdf 's for the four cosmological models for Zs = 1 and smoothing angles 
ranging from = 1' — 8'. Figures 9 and 10 show the dependence on Zs, by showing the pdf's for = 1' with Zs = 0.5 
and Zs = 2 respectively. Several qualitative features in P{k{9o)) can be understood from the results shown. As we decrease 
the smoothing angle, the peak of the pdf curves shifts towards negative values of k(9o), the peak height decreases and the 
distribution becomes more non-Gaussian. The reason is that when we decrease the smoothing angle we probe more non-linear 
scales in the underlying mass distribution and therefore the probability distribution for k(0o) is more non-linear. Increasing 
the source red-shift Zs has two distinct effects. It introduces more dark matter regions between the observer and the source 
which are virtually uncorrelated with each other. This has the effect of making probability distribution more Gaussian, even 
though the variance increases due to the increased path length. The peak height of the PDF decreases with increasing red-shift 
Za and P{k{9o)) becomes broader with the increase in the variance. 

For given smoothing angle and fixed redshift, the PDFs of the OCDM model are peaked most sharply, followed by that of 
the LCDM, TCDM and SCDM models. With the reduced convergence field ri{6o), P{ri{6o)) does not depend on the background 
geometry of the universe. We have noted that due to the presence of an extra |Kmin| in the denominator, the variance of the 
reduced convergence actually decreases with red-shift Zs in contrast to the variance of k,{9o). The skewness and other higher 
order moments for rj are also the same as for underlying mass distribution 5. It is therefore easy to understand that the PDF 
of the rj field is independent of background cosmological parameters and is very similar to the PDF of the density contrast 5. 
However the difference is that variance of 5 is much larger than unity, while that of rj is much smaller than unity, although 
the moments and hence the VPF associated with them are exactly the same. The cosmological dependence of enters 
through Kmin and {■rf{6o)). Since P(ry(&o))/|Kmi„| = P{k{9o)), the peak height of Pirj) is determined by the relative ordering 
of |Kmin|, for giveu red-shift. Thus we find that and {rfido)) play a very important role in the construction oi P{k{6q)). 

The approximation we have used to simplify the void probability distribution function of r] also gives us simple and new 
powerful method to compute the Sjv parameters for the convergence field, by relating them with with the Sn parameters of 
the underlying mass distribution as — Sn / {\nmin\)^~^ ■ 

The comparison of the analytical results with numerical simulations shown in figures 4-8 shows that there is a very 
good match, particularly for small smoothing angles. The exponential tail for positive values of k is well reproduced by our 
analytical results. The analytical PDF also reproduces correctly the exact position of the maxima and the sharp fall of the 
PDF for negative values of k. The relative ordering of the peaks for various cosmological models are also well reproduced. 
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Figure 9. Analytical predictions for the probability distribution function P(k), smoothed with angle 9o, as a function of k. The smoothing 
angle varies from 1' to 8'. The source red-shift is fixed at Zs = 0.5. The curves with increasing peak heights correspond to the SCDM, 
LCDM, OCDM and TCDM models. 



For large smoothing angles there is some disagreement between the analytical predictions and the ray tracing experiments. 
This is partially due to the fact that local initial spectral index for length scales making the dominant contribution for larger 
smoothing angles is different compared to that for the case of small smoothing angle. We have chosen uj — .3 for all our 
theoretical models, which correspond to initial power spectral index of n = —2. 

For larger smoothing angles length scales contributions to the one-point statistics of the convergence field come from 
quasi-linear length scales, hence the non-linear theory will not be accurate. In such cases one has to change the generating 
function G{t) to the perturbative generating function Q^'^ ij). With this change, we can extend all our results in a straight 
forward way. From the point of view of the Edgeworth expansion, the only change we need to make is in the values of ^jv 
parameters, where we have to replace the values obtained from the hyper-extended perturbation theory with ones obtained 
from the tree-level perturbation theory. However for very large smoothing angles [Oq > 1°) we have to keep in mind that the 
small angle approximation is a necessary ingredient in all our analytic calculations. 



5 DISCUSSION 

We have found an analytic expression for the PDF P[k{6o)) of the smoothed convergence field k{9q). This is a generalization 
of our earlier work on lower order cumulants. We found very good agreement of our analytic results with results from ray 
tracing experiments. We found that the PDF of the convergence field k(So) can distinguish different cosmological models. We 
also show that one can define a reduced convergence field rj{6o) which has the same reduced moments as its counterpart for 
the density contrast 5. 

While the variance of the convergence increases with source red-shift, the PDF itself tends to become more Gaussian. 
The differences in PDFs of different cosmological models become less prominent for very large source red-shift. From an 
observational view point this means that while the detection of weak lensing is easier from deep red-shift surveys, distinguishing 
cosmological models is easier when the source red-shift is smaller. We have not included the effect of noise due to intrinsic 




Figure 10. As in Figure 9, with source redshift Zs = 2. 



source ellipticity in our calculations. Several authors have focussed on weak gravitational lensing effects in the determination 
of cosmological parameters from SNela observations (e.g. Valageas 1999). Our results match with their finding for small 
smoothing angles. 

The study relating the PDF of the convergence field k(&o) with the PDF of the density contrast 5 can be easily extended 
to relate the Ke^ field with the bias associated with high peaks of the density contrast field 5. In recent studies (Munshi, Coles 
& Melott 1999a,b,c; Bernardeau & Scheaffer 1999) it has been shown that the generating function approach used here can be 
generalized to compute not only the bias associated with the over dense cells but can also provide valuable insight into the 
statistics of collapsed objects such as cumulants and cumulant correlators. Using the same analytical tools as we have used 
here, it can be shown that the statistical properties of high k(6'o) spots in convergence maps, e.g. cumualnts and cumulant 
correlators associated with them, can also be related to their counterpart for the density field 5. A detailed analysis of the 
signal to noise ratio is needed to determine the feasibility of estimating these quantities from observational data. 

The analytical expression for k used in this work is a line of sight integration of the density field 5 which relies on the Born 
approximation, i.e. it neglects higher order correction terms in the photon propagation equation. The error introduced by such 
an approximation in the quasi-linear regime has been studied by several authors (e.g. Schneider et al (1997). Perturbative 
calculations tend to show that this error is negligible for lower order cumulants, clearly such an analysis is not possible in 
the highly non linear regime. However since the tail of P{k{6o) contains information on all orders, a good match between 
theoretical predictions and the simulation results found by our study indicate that such corrections are negligible even in the 
highly non-linear regime. 
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